Non-singular solutions to Einstein-Klein-Gordon equations with a phantom scalar field 
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It is shown that the 4D Einstein-Klein-Gordon equations with a phantom scalar field (a scalar field 
with a negative sign in front of the kinetic energy term of its Lagrange density) has non-singular, 
spherically symmetry solutions. These solutions have a combination of features found in other 
O , spherically symmetric gravity plus field solutions. A stability analysis on these solutions indicates 

isJ '< they are unstable. 

U 
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I. INTRODUCTION 



ff^ ■ When one considers the system of gravity plus some field(s) (e.g. scalar field, gauge fields) there are a range of 
spherically symmetric solutions. When these extra fields are non-interacting, complex, scalar fields one finds the boson 
star solutions P, [3, Q which are prevented from collapsing because of the Heisenberg uncertainty principle. Allowing 
QQ the complex scalar field to have a self interaction of the X(p'^ type lead to boson stars with masses comparable to the 
. Chandrasekhar mass [1, in contrast to the much smaller mass of the non- interacting boson star solutions of 0, [1] ■ 
■ If one considers a real scalar field rather than a complex scalar then there are no non-singular, static solutions 
, [1, However it was also shown in that if one considers ghost fields (a real scalar field with a negative 

sign for the kinetic and potential parts of the scalar field Lagrangian density) one has non-singular, static solutions, 
but having a non-trivial topology. These ghost field solutions have a wormhole structure. This might have been 
expected since ghost fields may have a negative mass-energy density, and it is just such a form of matter that is 
needed in order to support wormholes fTT'|. A similar situation occurs when one considers wormhole-like solutions 
with a non- minimally coupled scalar field having a non-trivial topology (lo| . A final related system occurs with a 
phantom Born-Infeld field forming compact objects ( "gravastars" ) p^ . 

In this work we investigate spherically symmetric solutions to gravity plus a real scalar field. In contrast to previous 
work we consider a phantom field ^13j i.e. a real scalar field with a negative sign in front of only the kinetic energy 
part of the Lagrangian density. Phantom fields were originally proposed as an extreme form of dark energy which 
could explain the accelerated expansion of the Universe Phantom fields can violate the weak energy condition 
(WEC) and appear to be quantum mechanically unstable. (The WEC requires that the mass-energy density, e, and 
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pressure, p, of the field/fluid satisfy e + p > or w = p/e > —1). Any field or fluid with w < —1/3 will lead to 
gravitationally repulsion and is called dark energy. Phantom energy with w < — 1 is an extreme form of dark energy. 
A field or fluid with w > —1/3 leads to gravitational attraction. Nevertheless, as pointed out in 13] the supernova 
data seem to favor an extreme form of dark energy like phantom energy. The original proposal for phantom energy can 
be found in (see also ^ and references therein for related work in this direction) . There are also the so-called 

"quintom" models of dark energy [l^ with an equation of state that crosses the cosmological constant boundary and, 
correspondingly, violates the WEC. 

Phantom energy can lead to gravitational repulsion and therefore accelerated expansion of the Universe. This 
repulsion can provide a mechanism for the existence of static, regular solutions - the gravitational repulsion of the 
phantom field can be balanced by the gravitational attraction of normal matter. This is the physical mechanism 
behind the Bartnik-McKinnion solution [l^ to the Einstein- Yang-Mills system - the repulsion coming for the self- 
interaction of the Yang-Mills fields balances the usual gravitational attraction and stabilizes the configuration. As with 
the Bartnik-McKinnion solution the phantom field solution presented here is regular everywhere and asymptotically 
goes to Minkowski spacetime. Unlike the ghost field, kink solution of 0, with its non-trivial, wormhole topology, the 
present phantom field solution (like the Bartnik-McKinnion solution) represents a spherically symmetric, localized 
field configuration, and has a trivial topology. 

One unusual feature of the phantom field solutions examined here is that it has an overall negative mass. This is in 
contrast to the ghost field solution of 6] with its positive mass. Both phantom fields and ghost fields face difficulties 
when one tries to quantize them. Nevertheless the data on the accelerated expansion rate of the Universe favors 
something like phantom energy [l^ and so one should study the consequences of such unusual form of matter-energy. 
As far as we know phantom energy has been studied in a cosmological context where the phantom field varies only 
with time but does not vary spatially. In the present paper our phantom field varies spatially. Also strictly speaking 
our scalar field is only a phantom field (i.e. a field with w < —1) in some spatial regions, while in other regions it 
behaves as a gravitationally attractive field with w > —1/3. Physically this makes sense since a field which is phantom 
everywhere in space could not form some static spherically symmetric configuration but would disperse. 

There have been other works which have studied spherically symmetric solutions with phantom or ghost fields 
[IS HHj HI] ■ [13 black hole solutions in the presence of a phantom field were investigated. In ^24j the process of 
accretion of phantom energy onto a black hole was considered. In 25| a study was undertaken of spherically symmetric 
solutions for gravity plus some matter/fluid which violated weak energy condition. In this paper we replace the fluid 
by a scalar field. In regard to the negative mass of our solution there are other known negative mass solutions such as 
the negative mass black holes of jl^. However as with the ghost field kink solutions of Q these negative mass black 
hole solutions have non-trivial topology. 



II. EQUATIONS AND SOLUTIONS 



We chose our spherically symmetric metric to have the form 



(1) 



For the matter Lagrangian we have one real phantom scalar field with the Lagrangian 



L = 



R 



For the potential we take 



(2) 
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(3) 



where m is the mass of the scalar field and n is the coupling constant. This is the usual "Mexican hat" symmetry 
breaking potential. Note that in our terminology phantom field means a negative sign in front of the kinetic energy 
term only. By ghost field we mean a negative sign in front of both the kinetic and potential terms of the scalar field 
Lagrangian. With this definition it is phantom fields (not ghost fields) which can lead to w < —1 and accelerated 
expansion of the Universe. In what follows we set IQttG — 1. 
The energy-momentum tensor follows from ^ 



1 



(4) 
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The corresponding Einstein equations are 



(5) 



Note we are using the conventions of where the Ricci tensor is defined via the contraction of the first and third 
index of the Riemann tensor R^^ = R'^^v This is the reason for the overall + sign on the right hand side of ([5]) (see 
also (Hi). From ^ it follows 
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The wave equation for the scalar field Lp follows from variation of the Lagrangian ([2|) 



V — \ \ — e Lp' 



H + i(.'-A', 



+ 0. 

dip 



(7) 



In the above the ansatz functions are t and r dependent. Differentiation with respect to t (r) is denoted by a dot 
(prime). 

To investigate the static problem we drop the time dependence and express the metric functions v, A via new metric 
functions C(r),M(r) 



C{t) 



2M{r)' 




2M{r)' 


r 




r 



(8) 



(In the next section when we investigate the stability of our solutions we will again consider time dependent ansatz 
functions). In order to perform the numerical analysis we also introduce new dimensionless variables 



X = mr, M — 2mM, A — — - . 



Using these functions and variables the Einstein equations ([6]) become 
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and the equation for AI becomes 



M' = — x^ 
4 



Combining ^ and ^TU\i we find 



C 



Also the scalar field equation ([7]) takes the form 

1 / M' -M/x' 
^ 2 i 1-M/x 



M 



Thus we only needed to solve for M{x), ip{x) via pO)) (|12p since from pT|) C{x) was determined by ip{x). 
We solved equations (jlOp and (fT^ numerically using the following boundary conditions: 



(10) 



(11) 



(12) 



ip{0) = const, tp'{0) = 0, M(0) = 0. 



(13) 
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That is we demand regularity at the origin. Next choosing some A we found finite energy, non-singular solutions 
by numerically solving the system of equations ifTO]) and (|12p . For a given A this amounted to finding the correct 
initial value of the scalar field ip{0) which gave the desired well behaved solution. The problem reduces effectively to 
an eigenvalue problem since for a given A one had to search for the proper "eigenvalue" v?(0) in order to get a well 
behaved solution (regular, asymptotically flat) for the metric and scalar field ansatz functions. The flatness of the 
spacetime at infinity was provided by choosing the proper initial value for C(0) in equation (|lip i.e. we could chose 
C(0) such that C(oo) ^ 1. 

In figure [T] we show the graphs of the scalar field for two values of A. Generically the solutions for the scalar field 
looked similar - starting at some non-zero value at x = and decaying to zero as x — *■ oo i.e. (p{oo) — 0. For a normal 
scalar field this asymptotic value would be strange since for the "Mexican hat" potential this is a local maximum. 
Keeping in mind the unusual character of the phantom field it is not unexpected that the field should asymptotically 
go to a maximum rather than a minimum. The behavior of the metric functions e'''-^^ and e^^^-* for these values of A 
is shown in figure O Again the behavior of the metric functions is generally similar for different values of A - starting 
at some constant value at a; = they asymptotically go to 1 without ever becoming or negative. Thus the metric 
is non-singular and without horizons. 

The mass of the solutions was obtained from the asymptotic values of M{oo) from pH)) . The mass scale of the 
dimensionless quantity M{oo) is set by the mass m from the potential V{(p). The dependence of the ip{0) and M(oo) 
is shown in Table I and also in figure [3l From the figure and table one can see that the total mass goes to zero as 
A oo. Also ip{0) goes to zero as A — > oo. Af(oo) approaches zero from below while (p{Q) approaches zero from 
above. 



A. Asymptotic Behavior 



We now examine the asymptotic behavior of the solutions in order to get some analytic expressions for the solutions 
for large x and to compare to the numerical solutions. For this purpose, we write the asymptotic forms of the functions 
as: 

ifKidip, C « Coo + (5C, M « Moo + 6M, (14) 
where S(j),SC,6M <^ 1. Then the scalar field equation (fT^ becomes 

Sip" + -Sip' = Sip. (15) 

X 

We neglect terms of order {Sip)"^ , l/x^ and higher. In order to neglect the term K{Sip)'^ we need to additionally require 
that A is not too large. Although note that as A becomes large {Sip) goes more rapidly to zero. This last equation 
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TABLE I: The dependence of (/>(0) and M{oo) on A. 





A 


M(oo) 


7.1312051565 


0.1 


-23.8587 


5.09014118325 


0.3 


-3.68446 


4.35754 


0.5 


-1.85183 


3.478853 


1 


-0.821492 


2.697266 


2 


-0.391964 


2.288345 


3 


-0.256198 


1.3326834 


10 


-0.0755891 


0.6098228 


50 


-0.0152216 


0.4324704 


100 


-0.00756343 


0.1371201 


1000 


-0.000753043 





FIG. 3: The dependence of dimensionless mass M(oo) on A. pjQ 4. ^j^g equation of state w{x) when A = 0.3 (solid hne) 
Here small squares show calculated values of the mass and the j^jjfj = 10 (dashed hne) 
line is an approximation function with equation « const/ A^'^^ . 



has the asymptotic solution: 



S^^C,'^^. (16) 



Inserting this solution into (|10p gives the following equation and solution for SM: 

1„. , . . 1 



6M' ^~^clexpi-2x) ,5M « -C^ exp (-2a;). (17) 



Using these results in (fTT|) we have the following equation and solution for 6C: 

5C'«-ic^C^^^PlzM ^ SC^-lc^C'Ei{~2x), (18) 



where Ei is the exponential integral function which goes to from below as x — > 00. This shows that all perturbations 
{d(j), 6C, 6M) tend asymptotically to zero and that the ansatz functions as given in (|14p approach the values given by 
the numerical solution. 
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B. Behavior of the w parameter 

We now examine the behavior with respect to x of the effective equation of state w{x) = p{x)/e{x) of the scalar 
field (e(a;) and p{x) are respectively the energy density and pressure of the scalar field). From and ([5]) we have: 

T° = e{x) = -i _ + v{^), (19) 

n = -p{x) = 1 - + V{^). (20) 

The corresponding effective equation of state for the scalar field is: 

w{x) = —— = . 21) 

e{x) _i (i_M)^/2 + y(^) 

Using the numerical solution obtained earlier, figure |4] shows the behavior of w{x) for the case when A = 0.3 (solid 
line) and A = 10 (dashed line). From the figure one sees that there is some point x — x^ in which the denominator 
of (j2ip goes to zero (in the case under consideration this occurs around ~ 0.4). As x x^ from the left w — > —oo 
and from the right w +oo. Asymptotically w{x) tends to zero. In the range < a; < we have 

w < — 1 (w = — 1 at a; = 0) and the WEC is violated. Thus in this region our scalar field acts as phantom dark 
energy [29|. In the range a;* < a; < oo we have mostly gravitationally attractive matter with w > —1/3. For A = 0.3 
there is some region around x ~ 17 where w briefly dips down and just reaches w = —1. In this region the fleld is 
dark energy (but not phantom energy) and is gravitationally repulsive. However for both A = 0.3 and A = 10 the 
field is all or mostly gravitationally attractive. This provides a physical basis for the formation of these solutions: the 
gravitational repulsion coming from the regions with phantom energy {w < —1) or dark energy (—1 < w < —1/3) is 
balanced by the gravitational attraction coming from the regions where the field has w > —1/3. This is reminiscent 
of the Bartnik-McKinnion solution which owes its existence to the interplay with the repulsion of the Yang-Mills 
field against the attraction from gravity. In the present case the repulsion of the regions where the scalar field has 
w < —1/3 (and especially regions where w < —1) is balanced by the attraction of the regions where w > —1/3. 

There have been studies [30] where the equation of state of some fluid or field, as given by w, varies with time. In 
these studies there is some point in time where scalar field makes a transition and begins to violate the WEC and 
becomes a phantom field. In the present work our equation of state, as given by w, is spatially varying and it is this 
feature which gives rise to the static, spherically symmetric solutions. 

A similar result can occur even for a regular, real, scalar coupled to gravity. In [31] such a system was studied and 
it was found that there are regions where the weak energy condition is violated. This paper also found non-singular 
spherically symmetric solutions. 

III. STABILITY ANALYSIS 

We now study the dynamical stability of the above solutions against linear perturbations. We perturb the solutions 
of the system ^ by expanding the metric functions and scalar field function to first order as follows [11] [sl] 

v{t,x) ~ ^(^{x) + h'i{x) cos{ujt), 
X{t,x) ~ Xo{x) + Xi{x) cos{ujt), 
(p{t,x) — Lpo{x) + ip I (x) cos (tut)/ X. 

The index indicates the static background solutions of equations pUjl - lfT^ . (Expressions for i^oix), Xo{x) can be 
obtained from ([8]).) Then the first-order perturbation equations following from ([6]) and ([7]) are 

Ai = -l^Wi, (22) 



and the equation for tpi is 



V'l + 1(^0- K) f'l - Vo{x)vi + w^e^o-^'Vi = (24) 
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with the potential 



(1 - iK^l) 



(25) 



Introducing the new independent variable p 



dx 



we can rewrite equation (|24ll in a Schrodinger-like form 



where — e'^°^'^°Vo{x). If there is a negative eigenvalue cj^ then the solution will be unstable since then 

ipi ~ e*'^' will grow exponentially. To this end we examine the asymptotic behavior of the potential V^[a;(/9)]. Using 
the asymptotic expressions for the metric and scalar field ansatz functions from the previous section one finds 

V[x{p)] -e"'' = -1 as X oo. 

On the other hand, as x — s- we can expand the solutions into series in the following form 



Mo 



m^\p{^) [i-A^(o)2] x\ 




^(0) l-A^(O)^ 



2 4 

a; 



whence the expressions for i/q and Aq can be found (see equation dH)). Then, using ([SS]) . one finds that as x ^ 

l/[a;(p)] ~a + /3a;2, 

where 

a = 7 - C(0) [1 - 3A^(0)2 - 7] , /3 = C(0) {A^(0)2 [l - A^(0)2] - 7} , 

with 
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1 

12 



A 



<^(0)^ 



^(0)^ 



Depending on A and (^(0) the coefficients a and /3 can take values such that near x = one has a positive definite 
potential well or not. However since y[a;(p)] — > —1 as a; — > 00 one always has some negative uP' coming from continuum 
states. Because of these negative eigenvalues the solution is unstable. 

The presence of unstable solutions is not surprising. It is well-known [3^ that in 4D there are no stable soliton 
solutions to the Klein-Gordon equation with the "Mexican hat" potential (or any other type of potential). In Q 
it was found that adding gravity did not give stable solutions. Here we have added the additional ingredient of a 
phantom field but as the above stability analysis shows this also does not lead to stable solutions. 



IV. CONCLUSIONS 



By studying the system of 4D gravity plus a phantom scalar field (a real scalar field with a negative sign in front 
of the kinetic energy term in the Lagrangian density) we have found finite energy, regular solutions which have trivial 
topology. The solutions share some common features with previous solutions to systems of 4D gravity plus some 
field(s). Like the Bartnik-McKinnion solutions our solutions have finite energy, are regular everywhere, and have no 
horizons. The physical mechanism behind the present solutions is similar to that of the Bartnik-McKinnion solutions: 
the existence of the Bartnik-McKinnion solutions comes from the interplay of the repulsion of the Yang-Mills fields 
against the attraction of gravity. The present solutions arise from the balancing of regions where the scalar field a 
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phantom/dark equation of state (and is therefore gravitationally repulsive), against regions where the scalar field has 
the equation of state for ordinary matter-energy (and is therefore gravitationally attractive). 

The present solutions had some different features compared to the closely related ghost field, kink solutions of 
(for these solutions there was a negative sign in front of both the kinetic and potential terms in the scalar field 
Lagrangian density). These kink solutions had a non-trivial wormhole topology which resulted in their stability. In 
contrast the present solutions had a trivial topology, an everywhere regular metric and no horizons. However the 
linear stability analysis of the present solutions showed that they are unstable, as are the solutions of a regular scalar 
field plus 4D gravity and the higher Bartnik-McKinnion solutions (those with k > 2). One possibility toward finding 
stable phantom field solutions might be to alter the form of the potential V{(p). Note that asymptotically (p ^ 
which is an unstable equilibrium point for a normal scalar field. This behavior is due to the reversed character of 
a scalar field with a negative sign in front of kinetic energy term - it should go toward maxima of the potential. 
Unfortunately for the "Mexican hat" potential y = is only a local maximum; the global maxima is sd ip — ±oo. 
This suggests trying some potential which has a global maxima at some finite tp and with V(tp) finite at this point. 

The mass of our solutions turned out to be negative. For small A the magnitude of the mass was large and decreased 
to zero from below as A ^ oo (see figure [3] and table I). Given the unusual nature of the scalar field this could have 
been expected, and is similar to negative mass black hole solutions 26]. Despite the instability and negative mass 
features of the present solutions there is nevertheless a physical motivation for studying this system since experimental 
evidence jl3| favors the existence of something like phantom energy. 
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